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ABSTRACT 

In a laboratory experiment, water-ice aggregates are trapped in a vacuum chamber 
at a pressure of 2 mbar due to photophoresis and thermophoresis. The particles are 
located between a Peltier element at the bottom at 250 K and a reservoir of liquid 
nitrogen at the top at 77 K. Particle sizes vary between 20 /zm and a few hundred ^m. 
It is found that 95% of all the particles rotate about their vertical axis. A qualitative 
model is developed which explains why particles should mainly align to and rotate 
around the vertical. The results imply that rotation does not decrease the vertical 
strength of photophoretically driven motion in, e.g., protoplanetary discs. 

Key words: methods: laboratory - planets and satellites: formation - planet-star 
interactions - protoplanetary discs. 
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1 INTRODUCTION 

To verify and detail scenarios of protoplanetary disc evo- 
lution, of material transport and to interpret astronomical 
observations, it is crucial to understand the interaction of 
small dust and ice particles with their surroundings, i.e., gas, 
light, and other ice and dust particles. Therefore, it would be 
favourable if free particles were observed and manipulated 
in laboratory experiments. F ortunately, a t echni que to do 
so was recently developed bv iKelling^ et all (|201ll '). In their 
work, they show that micron-sized water-ice particles levi- 
tate under the Earth's gravity, trapped in stable positions 
by photophoresis and thermophoresis. Both phoretic forces 
are induced by a temperature gradient. In the first case, 
the temperature gradient exists across the particle (internal 
gradient), while the gaseous environment is at a constant 
temperature, in the latter case, the temperature gradient 
exists in the gaseo us environment (ext ernal gradient). In 
the experiments bv iKelling et al.l (|201ll ). photophoretic lift 
is induced by the thermal r adiation of a surface at 250 K. 

As already discussed bv lKrauss et al.l (|2007l ). levitating 
particles (in their case dust particles) tend to rotate. The 
same has been ob served for the levita ting ice particles in 
the experiment bv iKelling et al.l (|201ll '). Therefore, we here 
use the above mentioned trapping technique to perform a 
detailed study of the rotation of levitating ice particles which 
are under the influence of photophoresis, thermophoresis, 
and gravity. Rotation can have profound implications for 
the applicability of photophoresis in protoplanetary discs. 
The analysis of the rotation is important in this context as 
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it has been and still is frequently argued to be an obstacle 
to photophoresis. The argument is that photophoresis needs 
an internal temperature gradient which might be levelled 
out for a rapidly rotating particle causing the photophoretic 
strength to decrease. Therefore, the main aim of this paper is 
to observe particle rotation and to estimate the implications 
on the effect of photophoresis. 



2 EXPERIMENTAL SETUP AND FIRST 
OBSERVATIONS 

Ice particles are dropped onto a Peltier element cooled down 
to 250 K which is situated in a vacuum chamber at room 
temperature (see Fig. [T]). A reservoir of liquid nitrogen at 
77 K is placed above the Peltier element at a height of 
28 mm. The chamber is evacuated to roughly 2 mbar. As 
soon as the pressure decreases to below 10 mbar, single ice 
aggregates begin to rise from the surface of the snow layer. 
Their differen t traje ctories have been described in detail by 
iKelling et al.l (|201ll ). We here concentrate on the ice parti- 
cles which levitate halfway between the Peltier element and 
the reservoir of liquid nitrogen. We regularly observe the 
formation of groups of up to a few dozens of particles which 
are typically located above the edges of the Peltier element. 
While they are firmly trapped at a fixed height, they still 
rotate. We do not consider the crystal or aggregate structure 
of the ice particles. Instead, we simply classify the particles 
as water-ice. 

The particles vary in size between 20 ^ra up to a few 
hundred pm. This range is currently biased by the limited 
resolution of the optical system and also depends on the 
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Figure 1. Experimental setup: ice particles are trapped between 
a Peltier element at 250 K and a reservoir of liquid nitrogen at 
77 K due to photophoresis induced by the thermal radiation of 
the Peltier element. The chamber is evacuated to roughly 2mbar. 

particle sizes ejected from the snow placed on the Peltier 
element. 

To resolve the rotation, we obtain video sequences at 
502 frames per second. Figure [5] shows as an example one 
frame of the recorded video sequences (upper panel) and a 
rotation sequence of one of the aggregates (lower panel). 
About 95% of all the levitating particles are vertically 
aligned which means that they spin about their vertical axis. 
We do not find particles rotating around their horizontal 
axis. The remaining 5% of the particles show either no re- 
solvable rotation - they seem to be frozen halfway between 
the Peltier and the reservoir of liquid nitrogen - or they have 
an additional velocity component along the 2;-axis, i.e., they 
are moving up or down and generally show more complex 
trajectories. We note that in absolute numbers these are 
only a few particles. 

The particles rotating around their vertical axis can be 
divided into two classes: about half of the rotating particles 
have an additional component of motion as they move on a 
horizontal orbit, a rotation bound to the spin of the parti- 
cles. This is to be expected as the particles are typically of 
irregular shape (see model below). During a full orbit, they 
perform one turn around their vertical axis. Furthermore, we 
see many particles without the additional horizontal move- 
ment. However, this might be a resolution effect if the orbit 
is small. 

Our sample contains about 100 ice particles of differ- 
ent sizes and shapes. Within the resolution limit this in- 
cludes almost spherical particles as well as cylindric parti- 
cles, but mostly very irregular-shaped particles. Horizontally 
extended as well as vertically extended particles - both ro- 
tating around the vertical - occur (see Fig. [2]). 



3 ANALYSIS AND DISCUSSION 

We now discuss two parameters of the ice particles, i.e., their 
rotation frequency and their shape, in order to estimate their 
influence on the photophoretic strength. 




Figure 2. Upper panel: Levitating ice particles observed be- 
tween the Peltier element and the reservoir of liquid nitrogen at a 
pressure of l.Smbar. About 95% rotate about their vertical axis. 
Lower panel: Rotation of a single ice aggregate. The image sec- 
tion is 230 /.tmx 160 ^m. 



3.1 The dependence of photophoresis on the 
rotation 

In general, particle rotation can change the temperature gra- 
dients of illuminated particles. However, as long as the parti- 
cles rotate around the vertical (as they do here), always the 
same side faces the thermal radiation source at the bottom. 
Such a rotation does not change the temperature gradient 
in a coordinate system fixed to the particle (body fixed). 
The strength of photophoresis and thermophoresis along the 
direction of the source of light and the gas temperature gra- 
dient do then not depend on the rotation. Therefore, the 
fact that at least 95% of all the particles rotate around the 
vertical is a major finding here. It shows that rotation in- 
duced by phoretic forces is not important for the strength of 
photophoresis and thermophoresis along the vertical. This 
experimental result clearly shows that systematic rotation 
does not alter the capability of thermal radiation to lift par- 
ticles by phoretic forces. 

This might no longer be true if an additional side- 
ward illumination and therefore a horizontally directed pho- 
tophoretic force is added. In protoplanetary discs this might, 
e.g., be starlight illuminating particles at the surface of a 
fiared disc. The question then remains how the resulting 
photophoretic strength is influenced by the rotation around 
the vertical. In case of a strong photophoretic force the ro- 
tational state will certainly change, but this is subject to 
future work. If the vertical force is dominant, we might con- 
sider the simple case where rotation would still be tied to 
the vertical axis. In this case, the rotation frequency would 
be important. Fast rotation might level out temperature dif- 
ferences and reduce the photophoretic strength induced by 
the added source of light. On the other hand, slow rotation 
might lead to a transverse component of the force and ac- 
celerate or decelerate a particle, which leads to radial drift. 

To quantify such effects, the basic parameter is the ro- 
tation frequency. Therefore, we measured this parameter for 
the levitating particles. While we expect a wide range of fre- 
quencies due to the random nature of particle asymmetries, 
this range might still be conflned to and depend on the size 
of the particle. As a measure of the particle size we chose 
the horizontal, maximum extent of the ice aggregates be- 
cause the horizontal extent determines the moment of inertia 
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Figure 3. Large panel: the frequency of rotation of the ice par- 
ticles is plotted vs. their maximum horizontal extent, here called 
major diameter (both on logarithmic scales). Small symbols rep- 
resent all particles which show resolvable bound rotation. Large 
symbols represent all particles where a horizontal orbit does ei- 
ther not exist or is too small to be resolved. The uncertainties of 
both parameters are small (about the size of the large symbols) 
and for better visibility not plotted. Small panel: mean values 
of the frequency of rotation calculated for different ranges of the 
major diameter (10-20 /xm, 20-30 /jm, 100-260 /^m). 



around the vertical axis. The frequency plotted over the par- 
ticle size is shown in the large panel of Fig. |3] Small symbols 
represent all particles which show resolvable bound rotation. 
Large symbols represent all particles where a horizontal or- 
bit cannot be resolved. There is no clear trend visible. How- 
ever, the distributions of particles with and without notice- 
able bound rotation are ofTset with respect to their measured 
frequencies. While particles without bound rotation cover a 
range of frequencies from roughly 10 Hz to 100 Hz, particles 
with bound rotation need more time for one full turn and 
spin with a few up to 50 Hz. 

In the small panel of Fig. |3] we plot the mean of the fre- 
quency of rotation vs. the corresponding range of diameters 
(10-20/im, 20-30/im, 100 - 260 ^m). The plot suggests 
that large particles generally tend to have lower frequencies 
than smaller particles. 

We now compare the time a particle needs for a full rota- 
tion with the time that is need ed for conductive heat t rans- 
fer in the particle. According to lKrauss^ fc Wurml(|2005l ). the 
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time-scale is given by Thcat = , where p is the ice density, 
c is the heat capacity of the ice aggregates, and kth the ther- 
mal conductivity. For a typical particle radius of 10 fim with 
a density of 1000 kgm~^, a heat capacity of 1000 JkgK"^, 
and a thermal conductivity of about 1 W m~^ K~^, this char- 
acteristic time is 10~* s or 10000 Hz. This is sufficiently short 
to maintain a temperature gradient across the rotating ice 
particle. An experimental verification of this estimate is sub- 
ject to future work. 



3.2 The photophoretic alignment of a particle 

As shown by iKelling et al.l l|201ll ). four forces infiuence 
the behaviour of the ice aggregates: the thermo- and pho- 



tophoretic forces which accelerate the particles from the 
Peltier element towards the reservoir of liquid nitrogen, the 
gravitational force which opposes the phoretic forces, and 
gas drag. Gas drag caused by convection within the chamber 
is not important. However, gas drag is important as friction 
for particles moving upwards (which are not trapped) and ul- 
timately, gas drag also determines the rotation frequency as 
it provides damping. In equilibrium, the decelerating torque 
by gas drag has to balance the accelerating torque by pho- 
tophoresis and thermophoresis. We note that the difference 
in temperature inducing phoretic torques and the particle 
morphology inducing torques by gas drag have not neces- 
sarily to be around the same axis. However, within the res- 
olution of our experimental setup, 95% of all the particles 
are well described by a vertical rotation axis (see above). If 
the few particles behaving strangely were induced by such 
a misalignment, we would expect a more continuous dis- 
tribution, which we do not observe. One reason might be 
that damping is proportional to the velocity of a particle 
surface. Initially, only photophoretic torques act as acceler- 
ating force, which dominates the rotation. A non-aligned gas 
drag torque might therefore not dominate the motion of the 
particle. This is obviously the case in the experiments, but 
a theoretical treatment is beyond the scope of this paper. 

Thermophoresis and photophoresis are based on the 
same principle of interaction between gas molecules and dif- 
ferently heated particle surfaces. Therefore, the same physics 
of particle alignment should work. We then have a situation 
where gravity acts in one direction and a force depending 
on a temperature gradient directly opposes gravity. 

Rotation due to the different torques and forces is 
a problem which is highly complex and cannot easily be 
treated theoretically. Nevertheless, we here provide a simple 
view why particles should align to and rotate around the 
vertical (see Fig. |4]): in terms of basic physics, the motion 
of a particle can be divided into the motion of the centre 
of mass due to the total force and a rotation around the 
centre of mass due to a torque around the centre of mass. 
To understand particle rotation, we change its shape from 
a perfect sphere to a realistic, asymmetric particle in three 
steps. 

Step 1: often considered, but rarely realised by nature 
are homogeneous, spherical particles. Due to its symmetry, 
a homogeneous, spherical particle does not experience any 
torque and will not rotate (or at least, any initial rotation 
will quickly be damped by gas friction). Phoretic forces are 
directed upwards and for a levitating particle the phoretic 
forces balance gravity (Fig. 3] a). 

Step 2: as a next step on the way to a more realistic 
shape we consider a spherical particle where the centre of 
mass is not the geometrical centre of the sphere. Such a 
particle will also experience a phoretic force upwards which 
might be balancing gravity. However, in addition the parti- 
cle is subject to a torque (Fig. [4]b, left). This particle will 
rotate until the centre of the particle and the centre of mass 
are aligned to the vertical. In that configuration, the torque 
vanishes. The particle might oscillate a few times but due 
to gas damping, the oscillation will subside and the particle 
will be aligned eventually (Fig. |4]b, right). 

Step 3: adding surface structure to the particle leads to 
additional torques due to the fact that the phoretic forces 
are normal to the surface elements. Therefore, photophore- 
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initial conditions aiigned 




Figure 4. Rotation models of ice particles of different morphol- 
ogy. A temperature gradient is induced by a thermal source radi- 
ating from below. The centre of mass (CM) is indicated by a black 
dot, the temperature dependent photophoretic force is marked by 
slim arrows, the resulting forces by bold arrows. From top to bot- 
tom: (a) a homogeneous, spherical particle; (b) a spherical particle 
where the centre of mass is offset (left); this induces torque which 
causes rotation until the centre of the particle and the centre 
of mass are aligned (right); (c) a spherical particle with surface 
structure seen from the side (cl) and from above (c2); again, a 
state of equilibrium can be established (right). 



sis does no longer only act vertically, but the net effect also 
includes a horizontal component (Fig. |4]cl, left). Neverthe- 
less, also in this case a new configuration will be established 
where torques around the vertical vanish (Fig. [4]cl, right). 
However, a torque around a horizontal axis can prevail. If 
gravity is balanced by the vertical component of the pho- 
tophoretic force, a torque around the vertical due to irregu- 
larities of the surface will usually be left, for which we also 
indicate a top view here (Fig. |3]c2). 

As mentioned before, rotation around the vertical does 
not change the top and bottom parts of the particle. With 
respect to a coordinate system fixed to the particle, this 
torque remains constant as illustrated in Fig. |3] c and will 
accelerate the particle around the vertical axis. 

So far, we concentrated on the torque, but in addition to 
torques we also have to consider that non-sphericity induces 
a total force acting horizontally. If tuned in a way that the 
vertical component is balanced by gravity, only the sideward 
motion remains. If the particle is vertically aligned, this force 



becomes body fixed, i.e., is tied to the particle coordinate 
system. If the particle rotates, so will the direction of the 
horizontal force. This naturally leads to a bound rotation 
where the orbit of the centre of mass has the same frequency 
as the spin of the particle. 

Our model ignores any precession which might result 
from differences between angular momentum and rotation 
velocity and which might misalign the axis of the particle. 
Rotational damping by gas drag might also misalign the par- 
ticle as gas drag is not necessarily inducing a torque around 
the same axis. However, as long as the differences are small, 
these misalignment effects will be superimposed and domi- 
nated by the main phoretic torque. 

Slowly rotating bodies imply that the bound rotation is 
also slow and particle orbits can be larger. This is consistent 
with our observations that bound rotation (motion of the 
centre of mass) takes place preferentially for slow rotators. 
The fact that there is no visible sign of misalignment for 
bound rotators, indicates that gas drag by horizontal motion 
does at least to a certain extent not infiuence the alignment. 



4 CONCLUSIONS 

Our experiments show how particles rotate under the in- 
fluence of a photophoretic force, for which our model gives 
a plausible explanation. As the rotation axis aligns to the 
direction of illumination, the rotation, independent of its fre- 
quency, does not change the strength of photophoresis along 
the line of radiation. At the given (slow) rotation frequen- 
cies, temperature gradients appear to adjust rapidly enough 
to allow additional photophoretic forces from secondary ra- 
diation sources. However, a verification of this theory is sub- 
ject to future work. 

In general, our work supports the concept of pho- 
tophoretic transport of l evitating particles over the surface 
of protoplanetary discs ((Wurm fc Haack.2009a .b,'). 
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